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ABSTRACT 

We estimate the hkehhood of direct injection of supernova ejecta into proto- 
planetary disks using a model in which the number of stars with disks decreases 
linearly with time, and clusters expand linearly with time such that their sur- 
face density is independent of stellar number. The similarity of disk dissipation 
and main sequence lifetimes implies that the typical supernova progenitor is very 
massive, ~ 75 — 100 Mq. Such massive stars are found only in clusters with 
^ 10"^ members. Moreover, there is only a small region around a supernova 
within which disks can survive the blast yet be enriched to the level observed 
in the Solar System. These two factors limit the overall likelihood of supernova 
enrichment of a protoplanetary disk to ^1%. If the presence of short lived 
radionucleides in meteorites is to be explained in this way, however, the Solar 
System most likely formed in one of the largest clusters in the Galaxy, more than 
two orders of magnitude greater than Orion, where multiple supernovae impacted 
many disks in a short period of time. 

Subject headings: planetary systems: protoplanetary disks — planetary systems: 
formation — stars: formation 



1. Introduction 



Meteoritic evidence for live ^°Fe, a short lived radionucleide (SLR), in the early So- 
lar System strongly suggests that a supernova occurred shortly before the formation of the 
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planets (ITachibana et al.l 120061) . 30 years ago, and in the context of another SLR, ^^Al, 
Cameron &: TruranI (119771 ) proposed that supernova ejecta may be incorpo rated into plan- 
etesimals via the triggered collapse of the pre-solar nebula. More recently, iHester fc Desch 
( 120051 ) have suggested that SLRs are directly injected into newly formed protoplanetary 
disks. 

The evolution o f young stella.r clus ters and circumstellar disks has been well studied 
( iLada fc Lada I l2003l : iHaisch et al.ll200ll ). This allows us to quantify the likelihood that a 
supernova will occur close enough in time and space to a planet-forming disk to provide the 
inferred abundances of SLR in the Solar System. In this Letter, we address two questions: 
what is the likelihood of SLR enrichment of a disk by a supernova and what cluster properties 
maximize this likelihood? 

We summarize the relevant information on cluster and disk properties that we use in our 
model in §5J We find an analytic solution for the enrichment likelihood in a cluster in §3.11 
under the simplest assumption of a starburst, generalize to a finite formation period with 
a Monte Carlo simulation in §3.21 and consider the effect of multiple supernovae in ^ We 
calculate the overall enrichment probability and determine the cluster size that maximizes 
disk enrichment in ^ 



2. Parameters of the problem 



2.1. Cluster number distribution 



We adopt a cluster number distributi on, dNddN^ oc ^, that is consistent with 



both young, embe dded clusters ^ 3 Myr (ILada fc Lada I l2003l ). and older, optically vis- 
ible open clusters ( lElmegreen fc Efremovl 119971 ). The minimum cluster size is largely a 
matter of semantics for our calculations since small clusters are exceedingly unlikely to 
have supernovae within the maximum disk lifetime. The maximum cluster size, determined 
from the radio measurem ents of the ionizing luminosity of HIT regions, is A^*.max = 5 x 10^ 
jMcKee fc Williamsl ^9971 ) . 

There are proportionally more chances to find a star in a larger group so the differential 
probability that a star is found within a cluster with members, dP/dN^ oc l/N^,. We 
will plot probabilities as a function of the logarithm of the cluster size and define P/v, = 
dP/dln iV* Aln A^,, = PoAln iV*. The constant, Pq, is determined by the condition that the 
integral of Pn^ over al l clusters be equal to the total probability that a star is born in a 
cluster. iLada fc Lada I (120031 ) estimate that between 70 — 90% of all stars form in clusters 
and that 90% of these form in clusters with A^* > A^*min = 100. Taking an average of 80% 
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for the former gives Pq = 0.8 x 0.9/ ln(A^,,max/A^*,min) = 0.085. 



2.2. Cluster expansion 



Lada fc Lada I (120031 ) show that the number of detectable clusters declines with age and 
estimate that only about 10% of clusters survive as recognizable entities beyond 10 Myr. 
The spatial dispersion of nearby moving groups such as the TW Hya and f3 Pic associations 
show how quickly stars migrate away from their siblin gs. For the greater Scorpiu s-Centaurus 
region, with a spatial scale of 200 pc and age ~ 30 Myr (IZuckerman &: Song||2004j ) . the implied 
expansion speed is ~ 3 km s~^. 

All clusters at a given age, regardless of stellar number, have a similar avera ge stellar 
surface density, = N^/nrc^, where rr is the cluster r adius (lAdams et al. II2006I ). For the 
list of small clusters, A^^ ~ 30 — 1000, in iLada fc Lada I (120031 ) . the average surface density is 



100 stars pc ^ at an average age t ^ 3 Myr. The clusters in ICarpenterl (120001 ) . defined 



from 2M ASS data, have lower surface densities, T,^, ~ 30 stars pc~^. For massive star forming 
regions, iMcKee Sz Tan! (120031 ) find a much higher characteristic surface density, lgcm~^, 
corresponding to ~ 10^ stars pc~^ but in much younger objects, t <C IMyr. We assume a 



constant expansion speed such that 



Vexpt and S^, 



oc t ^, so the equivalent stellar surface 



density at 3 Myr is at least one, and possibly as many as two, orders of magnitude less. To 
bracket the possibilities, we consider the values, EsMyr = 30, 100, and 300 stars pc~^. 



2.3. Stellar mass function 



The stell ar initia l mass function (IMF) appears remarkably uniform over a range of 
cluster sizes ( iKroupal 1200 ll ). We use the Scalo mass function for which the number of 
stars with masses, M > 8Mq, that will become core collapse supernovae are a fraction 
fsN = 3 X 10"' ^ of all star s and follow a power law distribution, dN^/dM^ oc M~^^^°'\ with 
index a = 1.5 (IScaldll986l ). The cumulative distribution of supernova progenitors in a cluster 
containing A^^, stars is therefore 



N,{> M,) = fs^N, 



{MJM^Y - 1 



* {MjMiY - 1 

where the lower and upper limits of the progenitor distribution are taken to be = 8 Mq < 
< Mu = 150 Mq. The exist ence of an u pper limit is clear, although its actual value can 
only be statistically estimated (lFigerll2005l ). 
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2.4. Disk lifetimes and supernova timescales 



Mid-infrared surveys of stellar clusters, for which an average age can be determined, 
show that the fract ion of stars with dis ks, /d, is observed to decrease from unity at < 1 Myr to 
zero at td = 6 Myr (IHaisch et al.ll200ll ) . The decrease in disk fraction to zero is approximately 
linear with time, 

t>td. 



(2) 



UV radiation from O stars can ra pidly photoevaporate the outer radii of protoplanetary disks 
(e.g. IStorzer &: Hollenbachlll999l) but sufficient mass to form planet ary systems may remain 



bound to the star (jWilliams et al.l l2005l : lEisner &: Carpenterl 120061). Several of the clusters 



in the Haisch et al. survey and other similar studies (e.g. iMamajek et al.ll2004l ) contain O 
stars and the average disk fraction does not appear to be adversely affected, at least for 
luminosities comparable to Orion. Note also that whatever the star formation scenario, 
whether instantaneous, gradual, or induced, it is effectively incorporated into this empirical 
formalism. 



For supernova timescales, we use an empirical fit to the ISchaller et al.l (119921 ) stellar 
evolution models, log;^o ^sn = l-4/(log;^o ^sn)"^'^, where tsN is in Myr and Msn is in solar 
masses. 

If all stars in a cluster are coeval, the similarity of circumstellar disk and massive star 
lifetimes implies that no more than about half the disks remain when the first supernova 
occurs, even for the most massive progenitors. The least massive star that could explode 
within 6 Myr is Msn = 30 Mq, and is only likely to be found in clusters with A^* ^ 2600. If 
clusters are formed more gradually, some disks may exist at later times and slightly lower 
mass progenitors may play a role (§3.2). 



2.5. Proximity to supernova blast 



Disk enrichment places tight constraints on spatial scales too. A minimum mass solar 
nebula (MMSN) disk with r adius 100 AU ar ound a solar mass star will be stripped by su- 
pernova ejecta within 0.2 pc (IChevalierl |2000| ) . Matching the abundances of the ejecta with 
the meteoritic recor d, however, requires that disks lie w ithin 0.22 pc for a 25 progenitor 
( iLoonev et al.l 120061 ) and 0.3 pc for a 40 Mq progenitor (lOuellette et al.ll2005l ). We find that 
even larger masses are more likely sources of enrichment due to their shorter main sequence 
lifetimes. These can enrich a larger volume, out to ~ 0.4 pc for a 100 M© progenitor. We 
therefore consider a radial range, < r < rs + r^,, where = = 0.2 pc within which disks 
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can survive the supernova and be enriched to the level observed in the Solar System. The 
number of stars in this "enrichment zone" depends on the cluster density profile and size. 

^Observations of Orion ( Hillenbrand fc Hartmannlll998 ) and other clusters (e.g. Muench et al, 



20031 ) show that the number of stars per unit area declines approximately inversely with an- 
gular distance from the center (albeit with some significant substructure in some cases). The 
inferred stellar volume density profiles are therefore approximately inverse square, n^: oc 
and the number of stars increases linearly with radius, N^{< r) = N^{r/r^, where is the 
cluster radius. The most massive stars in a cluster are generally found near its center 
( iLada fc Lada I l2003l ). Assuming that this is the case for the supernova progenitor, the 
fraction of stars in the enrichment zone is 



fe{t) 



N^{< rs + re) - iV,(< Ts 



(3) 



The stellar motions in a cluster have characteristic value, v, 



exp 



rjt = (iV,/97rS3Myr)'/'. 

For example, a typical cluster in iLada fc Lada I (l2003l ) with N^, = 10^, SsMyr = 100 pc~^, 
has fexp = 0.6 km s~^ but the velocities are higher in larger clusters with similar surface 
densities. To estimate the fraction of stars within the enrichment zone, we assume that 
the cluster maintains its inverse square density profile, and, when considering the effect of 
multiple supernovae (§!]), that stars move independently through this zone. 



3. Calculation of the enrichment likelihood 



3.1. Cluster formation in a starburst 



Under the assumption that all the stars in a cluster form at the same time, i.e. in a 
starburst, the first supernova will be the most massive star and the enrichment likelihood 
can be calculated analytically. 

The conditional probability that the most massive star in a cluster has mass Msn is 
the expected number of stars of this mass times the probability that there are none more 
massive (see 



Williams fc McKed 119971 ) 



[{MjMiY - 1]M, V^sN 



l+a 



-M*{>Msn) 



(4) 



If we further assume that the disk lifetime is independent of cluster location then the 
fraction of disks that exist within the enrichment zone can be separately factored as /d/e- 
Evaluating this at the time of the supernova and integrating over all possible progenitors 
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gives the likelihood that a disk in a cluster of a given size is enriched with SLR at the level 
observed in our Solar System, 

r-Mu 

-PSLR|N, = / /d(^SN) fe{tm) PusnW* dMstq. (5) 
JMi 

Figure [H shows that Pslr|n, has a similar form independent of EsMyr with a broad 
maximum centered on N^, ~ 8000. Small clusters are unlikely to have supernovae within 
the disk lifetime and large clusters expand more rapidly so relatively few stars lie within 
the enrichment zone. The absolute likelihood increases for higher surface densities since 
the expansion speed will be lower, and the enrichment fraction higher, for a given cluster 
number. 



3.2. Cluster formation over a finite duration 

The starburst assumption is a simplification: young clusters contain protostars in a 
range of evolutionar y states suggesting that they are built up over a period tgf ~ 1 Myr 



flLada fc Lada 1120031 ) 



For a prescribed star formation rate (SFR), N^,{t), the average disk fraction in the cluster 
is /d(t) = /q /d(t — t') N*{t') dt' /N:^. The time dependence of the SFR in clusters is unknown. 
The simplest assumption is that it is constant, = A^*/tsf- In this case, for t > tgf, 

1 - (t - tsf/2)Ad t < td, 
fdit) = { (tsf + U- tf/2UU U<t<U + tsf (6) 
t>U + tsf. 

We proceed by randomly sampling a star from the IMF at each time interval At = 
The birth time is added to the main sequence lifetime for each star with M^, > 8 Mq to 
determine when the first supernova occurs. At this time, ti, Pslr|n* is calculated as the 
product of the disk fraction above times the enrichment fraction, where the cluster size 
rc = "^exp^i = (A^*/97rE3Myr)"'^^^^i pc. As ti is a stochastic variable, we average over 10^ 
simulations for each cluster number. 

Figure [T] shows that -Pslr|n* is relatively insensitive to tgf , and much more dependent 
on SsMyr- Peak probabilities are slightly lower for a finite formation time as the higher disk 
fraction at the time of the supernova is offset by the larger cluster size and lower enrichment 
fraction The extended period over which a supernova can impact circumstellar disks allows 
for lower mass progenitors, however. The average supernova mass is 74 Mq for tsf = 1 Myr 
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compared to 98 Mq for the starburst. In each case, such high masses are favored because of 
the short disk hfetimes and are only hkely to be found in very large clusters, A^.,, ^ 10^. 



4. Multiple supernovae 

Large clusters, A^* ^ 5000, should have more than one supernovae within the maximum 
disk lifetime, + tsf- We model the effect of additional supernovae by summing enrichment 
probabilities using the same Monte Carlo formulation as in §3.21 The assumption here is 
that cluster dynamics move new star-disk systems into the enrichment zone by the time of 
the next supernova. Two relatively small corrections are made: first, the fraction of possible 
disks that may be enriched is decreased by the fraction that are stripped by the previous 
supernova, Vs/rc] second, the fraction of disks in the enrichment zone, /e, is decreased by 
Vexp^ts^/^e if this is less than 1, where AtsN is the interval between supernovae, to allow 
for migration of new disks into this region. In practice, only at most a few percent of disks 
are close enough to a supernova to be destroyed and, except for the very largest clusters, 
the interval between supernovae is large enough that the enrichment zone is continually 
refreshed. 

The results are shown in Figure [2] both for starburst and extended SFR scenarios, 
ts{ = 0,3 Myr. The enrichment likelihood continually increases to the largest clusters where 
many tens of supernovae can contribute. The inset shows the total number of supernovae 
within td + tgf and the number that enrich half of the overall total for a fiducial tgf = 1 Myr. 



5. Implications 

The overall probability of supernova enrichment for any star is determined by integrat- 
ing the enrichment likelihood, or conditional probability of enrichment given cluster size, 
over the cluster number distribution, Pslr = -Pslrin* -Po c^ln A^*. This "Galactic 

enrichment likelihood" is plotted versus a wider range of average surface densities in Fig- 
ure [31 The incorporation of the cluster number distribution adds some uncertainty to the 
absolute numbers but the general form, in particular the strong dependency on T^^uyr and 
near-independence on tgf, is inherited from Pslr|n, and is robust. As the surface density 
decreases with time, the cluster will become harder to identify. A rough estimate of the 
cluster lifetime, tc, for which S(tc) = 3starspc~^, comparable to the field star density, is 
shown on the upper axis. Given that ~ 90% of all clusters do n ot survive beyond 10 Myr, and 



that planetary systems are disrupted in longer lived systems (lAdams fc Laughlirul200ll ). we 



conclude that supernova enrichment of protostellar disks is a highly unlikely event, affecting 
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< 



1% of all stars in the Galaxy. 

How, then, to explain the presence of ^°Fe in the early Solar System? If, however un- 
likely, it was injected into the protoplanetary disk from a nearby supernova, then a Bayesian 
estimate of the most likely cluster size is -Pn.|slr = -Pslr|n,-P/v,/-Pslr- Because P/v* is in- 
dependent of A^* and Pslr is a normalization factor, the conditional probability of stellar 
number given enrichment is directly proportional to the enrichment likelihood plotted in 
Figures [H [2l This strongly favors our Solar System' s origin in the largest clusters in the 
Galaxy, exemplified by NGC3603 (IMoffat et al.lll994l ). more than two orders of magnitude 
more luminous than Orion, where multiple supernovae can potentially enrich many disks. 

A key assumption in this conclus i on is that the disk fraction is independent of cluster 
location. Recent work by iBalog et al.l (120071 ). however, shows that /d is a factor 2 — 3 lower 
than equation[2] in the central 0.5 pc of NGC2244, a cluster intermediate in luminosity be- 
tween Orion and NGC 3603. This would decrease -Pslr|n* by the same amount, lessen or even 
nullify the increase toward large in Figure [2], and only strengthen the conclusion that the 
direct injection of supernova ejecta into a protoplanetary disk is a very unlikely event. More- 
over, if one or more massive stars is the source of other SLRs, particularly ^^Al, in the early 
Solar System, then their implied scarcity in other planetary systems r nay have important 
impli cations for the thermal history of planetesimals in those systems (IGhosh fc McSween 
1998h . 



This work was motivated by engaging converstaions with Sasha Krot and Gary Huss 
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REFERENCES 

Adams, F. C., & Laughlin, G. 2001, Icarus, 150, 151 

Adams, F. C., Proskow, E. M., Fatuzzo, M., & Myers, P. C. 2006, ApJ, 641, 504 

Balog, Z., MuzeroUe, J., Rieke, G. H., Su, K. Y. L., & Young, E. T. 2007, ApJ, in press 
( [astro-ph/070174l [) 

Cameron, A. G. W., & Truran, J. W. 1977, Icarus, 30, 447 

Carpenter, J. M. 2000, AJ, 120, 3139 

Chevalier, R. A. 2000, ApJ, 538, L151 



-9- 



Elmegreen, B. G., & Efremov, Y. N. 1997, ApJ, 480, 235 

Eisner, J. A., & Carpenter, J. M. 2006, ApJ, 641, 1162 

Figer, D. F. 2005, Nature, 434, 192 

Ghosh, A., & McSween, H. Y. 1998, Icarus, 134, 187 

Haisch, K. E., Jr., Lada, E. A., & Lada, C. J. 2001, ApJ, 553, L153 

Hillenbrand, L. A., & Hartmann, L. W. 1998, ApJ, 492, 540 

Hester, J. J., & Desch, S. J. 2005, Chondrites and the Protoplanetary Disk, 341, 107 

Kroupa, P. 2001, MNRAS, 322, 231 

Lada, C. J., & Lada, E. A. 2003, ARA&A, 41, 57 

Looney, L. W., Tobin, J. J., & Fields, B. D. 2006, ApJ, 652, 1755 

Mamajek, E. E., Meyer, M. R., Hinz, R M., Hoffmann, W. F., Cohen, M., & Hora, J. L. 
2004, ApJ, 612, 496 

McKee, C. F., & Tan, J. C. 2003, ApJ, 585, 850 

McKee, C. F., & Williams, J. P. 1997, ApJ, 476, 144 

Moffat, A. F. J., Drissen, L., & Shara, M. M. 1994, ApJ, 436, 183 

Muench, A. A., et al. 2003, AJ, 125, 2029 

Ouellette, N., Desch, S. J., Hester, J. J., & Leshin, L. A. 2005, Chondrites and the Proto- 
planetary Disk, 341, 527 

Scalo, J. M. 1986, Fundamentals of Cosmic Physics, 11, 1 

Schaller, G., Schaerer, D., Meynet, G., & Maeder, A. 1992, A&AS, 96, 269 

Storzer, H., & HoUenbach, D. 1999, ApJ, 515, 669 

Tachibana, S., Huss, G. R., Kita, N. T., Shimoda, G., & Morishita, Y. 2006, ApJ, 639, L87 
WiUiams, J. P., Andrews, S. M., & Wilner, D. J. 2005, ApJ, 634, 495 
WiUiams, J. P., & McKee, C. F. 1997, ApJ, 476, 166 
Zuckerman, B., & Song, I. 2004, ARA&A, 42, 685 



-10- 



This preprint was prepared with the A AS macros v5.2. 



- 11 - 




Fig. 1. — The enrichment hkehhood from a single supernova event versus stellar number for 
three values of the average surface density at 3 Myr. The shaded regions show the variation 
with cluster formation time from tgi — — 3 Myr, with the starburst scenario indicated by 
the thicker line. 
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Fig. 2. — The enrichment hkehhood as a function of stellar number, as in Figure [T|, but for 
the case of multiple supernovae. The inset shows the total number of supernovae (solid line) 
per cluster and the number at which half the final enrichment fraction was reached (dashed 
line) for the case tgf = 1 Myr. 
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Fig. 3. — The dependence of the overall likehhood of SLR enrichment, summed over all 
cluster sizes and allowing for stars not formed in clusters, as a function of the average stellar 
surface density at 3 Myr. The variation with cluster formation time from tgf = — 3 Myr is 
indicated by the shaded region with the starburst scenario shown as the thicker line. 



